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We present numerical calculations of lamellar phases of block copolymers at patterned surfaces. 



We model symmetric di-block copolymer films forming lamellar phases and the effect of geomet- 

O ■ 

rical and chemical surface patterning on the alignment and orientation of lamellar phases. The 



calculations are done within self-consistent field theory (SCFT), where the semi- implicit relaxation 



scheme is used to solve the diffusion equation. Two specific set-ups, motivated by recent exper- 
iments, are investigated. In the first, the film is placed on top of a surface imprinted with long 



chemical stripes. The stripes interact more favorably with one of the two blocks and induce a 
perpendicular orientation in a large range of system parameters. However, the system is found to 
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iy\ | be sensitive to its initial conditions, and sometimes gets trapped into a metastable mixed state 

composed of domains in parallel and perpendicular orientations. In a second set-up, we study the 



film structure and orientation when it is pressed against a hard grooved mold. The mold surface 
prefers one of the two components and this set-up is found to be superior for inducing a perfect 
perpendicular lamellar orientation for a wide range of system parameters. 



I. INTRODUCTION 

Block copolymers (BCP) have been studied extensively in the last few decades due to their 
special self-assembly properties giving rise to interesting mesophases in the sub-micrometer 
to nanometer range, as well as to their numerous applications, where desired properties can 
be tailored by specific chain architecture [1-8]. 

Bulk properties of BCP are rather well understood and, in recent years, much effort was 
devoted to understand thin films of BCP. One potential application is the use of thin films 
of di-block copolymers as templates and scaffolds for the fabrication of arrays of nanoscale 
domains, with high control over their long-range ordering, and with the hope that this tech- 
nique can be useful in future micro- and nano-electronic applications. Recent experiments 
include using chemically [9-18] and physically [19-22] patterned surfaces, which have pref- 
erential local wetting properties for one of the two polymer blocks. The orientation and 
alignment of lamellar and hexagonal phases of BCP were investigated, and, in particular, 
their transition between parallel ('lying down') and perpendicular ('standing up') orienta- 
tions. Another useful method is the use of electric fields to orient anisotropic phases of BCP, 
such as lamellar and hexagonal, in a direction perpendicular to the solid surface [23-31]. 

In this paper, we present self-consistent field theory (SCFT) calculations inspired by 
recent experiments on patterned surfaces [10-14, 32-34]. Our main aim is to analyze what 
are the thermodynamical conditions that facilitate the perpendicular orientation of BCP 
lamellae with respect to the underlying solid surface, and how the lamellar ordering can be 
optimized. Two specific solid patterns and templates are modeled. The first is a planar solid 
surface that has a periodic arrangement of long and parallel stripes preferring one of the two 
blocks, but otherwise is neutral to the two blocks in its inter-stripe regions. We show that 
this experimentally realized surface pattern [10-14] enhances the perpendicular lamellar 
orientation. The second surface pattern is motivated by recent Nanolmprint lithography 
(NIL) experiments [14, 32-34]. This is a high-throughput low-cost process which has the 
potential of reducing the need for costly surface preparation. Here, a hard grooved mold is 
pressed onto a thin BCP film at temperatures above the film glass-transition and induces 
perpendicularly oriented lamellae. Within our model we show that, indeed, the grooved 
surface does enhance the perpendicular orientation of lamellae. 

The SCFT model that we use in the BCP calculations has several known limitations. It 



is a coarse-grained model and, as such, can only describe spatial variations that are equal 
or larger than the monomer size (the Kuhn length). Our calculations provide the thermo- 
dynamical equilibrium, or local minima of the film free-energy in presence of geometrical 
constraints. Therefore, important structural details induced by hydrodynamic flow and film 
rheology as occurring during sample preparation are not described by the model. 

In the present work we limit ourselves to three-dimensional systems that are translation- 
ally invariant along one spatial direction. This is applicable when the BCP film is put in 
contact with surfaces having long unidirectional stripes or grooves. Extensions of the present 
work to more complex three-dimensional systems with two-dimensional surface patterns will 
be addressed separately in a follow-up publication. 

The outline of our paper is as follows. In the next section two system set-ups are intro- 
duced; a chemically striped surface and a grooved mold and their effect on orienting lamellar 
BCP films is presented. In Sec. Ill we describe our SCFT model and how its equations are 
solved numerically. In Sec. IV our results are presented for the two types of experimental 
set-ups. Finally, in the last section we discuss the model predictions and their connection 
with experimental findings. 

II. THE BCP FILM DESIGN 

We consider a melt of A-B di-block copolymer (BCP) chains composed of n chains, each 
having a length N = N^+Nb in terms of the Kuhn length a that is assumed, for simplicity, to 
be the same for the A and B monomers. Hence, the A-monomer molar fraction / = N^/N 
is equal to its volume fraction. In addition, hereafter we concentrate on symmetric di- 
BCP, Na = Nb having / = 0.5. The symmetric BCP yields thermodynamically stable 
lamellar phases of periodicity £q, as the temperature is lowered below the order- disorder 
temperature (ODT) [35]. At shallow temperature quenches, simple scaling arguments [36] 
used in the weak segregation limit show that the lamellar period £q is proportional to R g , 



the chain radius of gyration, £ ~ R g = JNa 2 /6 ~ N l l 2 . For deep temperature quenches 
well below the ODT, the strong segregation theory [37] yields more stretched chains as 
£ ~ N 2 / 3 > Rg ~ N 1 ' 2 . 

The BCP film has total volume Q and lateral area A, so that its thickness is L = Q/A. 
In some experimental set-ups the BCP film is bounded by two planar solid surfaces, and its 



thickness L is a constant. In other set-ups [10-14] the film is spin coated on a solid surface 
with a free polymer/air interface on its top, so that the thickness can vary spatially. In yet 
another set-up used in Nanolmprint lithography (NIL) experiments [14, 32-34], a grooved 
mold is pressed against the film and the film penetrates into the mold. As the film profiles 
inside the mold varies considerably in height, L is only the film average thickness. 

We will consider only surface features along one spatial direction (chosen to be the x- 
direction), and assume that the system is translationally invariant along the second surface 
direction (the y-direction). Hence, the film volume Q (per unit length) has units of length 
square, while the surface area A has units of length. The third spatial direction, the z one, 
is taken to be perpendicular to the surfaces. This allows us to carry out the numerical calcu- 
lations only in the (x,z) two-dimensional plane, and represents a considerable simplification 
from the numerical point of view. 

The situation where a thin BCP film is placed in contact with a flat and uniform surface 
(or is sandwiched between two flat surfaces) was modeled by several authors [36, 38-46]. 
Two main features are apparent when the film behavior is compared to that of bulk BCP. 
The first effect is the film confinement. When L differs from the natural periodicity £q, the 
chains need to be stretched or compressed as the film is incompressible and space filling. 
The film free-energy shown in Fig. 1 is a function of the thickness L, and is obtained within 
our SCFT scheme (see below), and agrees well with previous results [38, 43]. 

The main effect of the confinement between the two bounding surfaces is the existence 
of free-energy minima at integer or half-integer values of L/£o = |, 1, |, ... corresponding to 
film thicknesses where we can fill an integer or half-integer numbers of A-B parallel layers 
in between the two surfaces. The overall trend for the film free-energy is to converge toward 
the bulk value F\, as: F — Fj, ~ 1/L. 

The second feature is the possibility to induce a parallel to perpendicular transition of 
the lamellae by changing the strength of the surface preference, Au. This can be seen in 
Fig. 2 where the parallel to perpendicular phase diagram is plotted in the Aw - L plane, 
within our SCFT scheme. When a strong surface preference towards one of the two blocks 
is included, the lamellae tend to orient in a parallel direction, while for neutral (indifferent) 
surfaces or weak preferences, the perpendicular orientation is preferred as the lamellae can 
assume their natural periodicity £o in this orientation for any thickness L. Note also that 
the transition occurs at Au = for integers or half-integers values: L/£ = |, 1, |, ... as was 
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FIG. 1: The difference of the dimensionless film free-energy from its bulk value, (F — Fb)L/(nkBT£o), 
as function of the film rescaled thickness, L/£q, where £q is the lamellar periodicity, k B T is the thermal 
energy and n the total number of chains. The lamellae in the film are parallel to the two flat bounding 
surfaces (the L|| state). The surface preference is Au = 1 for the bottom surface and Au = for the 
top, and Nx = 20. 

argued above. These results agree well with those reported in Refs. [36, 38, 43]. 

In the remaining of the paper we will address in detail the question of how it is possible 
to better control the relative stability of parallel and perpendicular phases of lamellar BCP 
films. And, in particular, how can the stability of the perpendicular phase be increased for 
a larger range of film thicknesses and surface characteristics. 

III. THEORETICAL FRAMEWORK 



Since the system is translationally invariant in the y-direction, we treat it as an effective 
two-dimensional system. The free energy for such a di-block copolymer (BCP) film confined 
between the two surfaces is 
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FIG. 2: Phase diagram for intermediately segregated (Nx = 20) symmetric di-block lamellar phase in 
terms of the film thickness L/£o vs the surface field preference between the two blocks, Am = ua — ub- 
The grey area indicates the perpendicular lamellar region (L^) and the white region above it is the 
parallel state (Lii)- The film is confined in the z-direction between two parallel and flat surfaces. The 
bottom surface uniformly attracts one of blocks (Am > 0), while the top one mimics the free (and 
neutral) surface, Au = 0. 



where each of the n BCP chains is composed of N = A^ + A# Kuhn segments of length 
a, and the Flory-Huggins parameter is x- The dimensionless volume fractions of the two 
components are defined as <pA( r ) = (J>a(x,z) and 0#(r) = <j>b(x,z), respectively, whereas 
UJ i{ r )i J = A, B, are the auxiliary fields coupled with <ftj(r), and Qc is the single-chain 
partition function in the presence of the ua and ujb fields [see Eqs. (4)-(5) below for more 
details]. The third term represents a surface energy preference, where ua and «g are the 
short-range interaction parameters of the surface with the A and B monomers, respectively. 
Formally, ua{t) and ub{t) are surface fields and get non-zero values only on the surface(s). 
Finally, the last term includes a Lagrange multiplier i](r) introduced to ensure the in- 
compressibility condition of the BCP melt: 



4>a{t) + 4>b{t) = 1 for all r e VL 



(2) 



By inserting this condition, Eq. (2), in the surface free energy of Eq. (1), the integrand 
becomes ua4>a + ub4>b = (ua — ub)<Pa + Ub- Hence, Au(r) = ua(j : ) — u B (r) is the only 
needed surface preference field that will be employed throughout the paper. 
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Using the saddle-point approximation, we obtain a set of self-consistent equations 
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where the incompressibility condition, Eq. (2), is obeyed, and the single-chain free energy 
Qc is: 

Qc = ^Jd 2 rq A (r,N A ) (4) 

The two types of propagators qj (r, s) and qj (r, s) (with j = A, B) are solutions of the 
modified diffusion equation 

^^ = J V 2 ^ (r, 8 ) - ^(r)g,(r, s) (5) 

with the initial condition q A (r, s—0) —q B {r, s—0) = 1, q A (r : s—0) — q B (r : N B ) and 
q B (r, s—0) = q A (r, N A ), where s is a conveniently defined curvilinear coordinate along 
the chain contour. This diffusion equation is solved using reflecting boundary conditions 
at the two confining surfaces (z = and z = L): dq/dr\ z=Q = and dq/dr\ z= L = 0, while 
periodic boundary conditions are used in the perpendicular direction. 

Hereafter, we rescale all lengths by the natural periodicity of the BCP, £o — 4.05R g , [47] 
where R g is the chain radius of gyration R 2 = Na 2 /6. Similarly, s is rescaled by N, yielding 
r — >■ r/£ , s — » s/N, x —> Nx, Wj(r) —¥ Nu>j(r) and Uj(r) — > Nuj(r) with j =A or B. With 
this rescaling, we rewrite the self-consistent equations as: 

wa(t) = X<t>B{r) - u A (r) + r}(r) (6) 

wbW = X<pA{r) - u B (r) + rj(r) (7) 

1 rf + 

4>A{r) = 7— / ds q A (r, s)q A (r, f - s) (8) 
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^^ = (^) 2 V 2 g,(r, S )-c,(r)g,(r, S ) (10) 



where / = Na/N, Oi = £1/£q and Qc = (l/^i) I d 2 r q A (r, f). Note that the incompressibil- 
ity condition, Eq. (2), together with Eqs. (6) and (7) can be used to obtain the Lagrange 
multiplier i](r) 

V( r ) = ^(ua + ojb-X + ua + ub) (11) 

With the rescaled variables, we define now a rescaled free energy: 

1 



\nQc - ^- / d 2 r [wA(r)0A(r) + u B (r)^ fl (r)] 
+ — |d 2 r?7(r)[0 A (r) + <Mr)-l] (12) 

The above self-consistent equations can be solved numerically in the following way. First, 
we guess an initial set of values for the auxiliary fields u>j(r). Then, through the diffusion 
equations, Eq. (10), we calculate the propagators, qj and qj. Next, we calculate the monomer 
volume fractions <pj from Eqs. (8)-(9) and the Lagrange multiplier r)(r) from Eq. (10). We 
can now proceed with a new set of values for 0Oj(r) obtained through Eqs. (6)- (7), and 
this procedure can be iterated until convergence is obtained by some conventional criterion 
described below. 

We use the semi-implicit relaxation scheme [48] to solve the two-dimensional modified 
diffusion equations, Eq. (10). Our convergence criterion is based on the incompressibility 
condition. For perfect structures such as parallel or perpendicular lamellae, the maximum 
allowed deviation between the sum of the A and B densities and unity, |1 — 0^(r) — 0s(r)|, 
is 10~ 4 , whereas for the mixed Lm phase (see below), it is around 10~ 2 . As mentioned 
above we rescale all lengths by the natural periodicity of the BCP, £q — 4.05R g , and the 
curvilinear coordinate, s, by the total number of monomers in one chain, N. The spatial 
discretization in the x-direction is 0.05 (in units of £q), while in the z-direction it is 0.025. 
The discretization of the s variable is 0.02. For all the presented results, the free energy 
changes in the last few iteration steps are less than 10~ 4 in units of /c^T/chain after the 
first 1,000 iterations and decreases to 10 -6 after additional 4,000 iterations. Note that since 
we work at a mean-field level (SCFT), it would not be of advantage to further refine the 
convergence of the free energies to a higher accuracy, since we neglect anyway quadratic 
fluctuations that might give larger corrections. 



IV. RESULTS 

We present now the numerical results for symmetric di-block films (/ = 1/2) at various 
patterned surfaces. The natural periodicity of the BCP, £q, is chosen for all the numerical 
calculations to be 50 nm. This value roughly corresponds to values used in several experi- 
mental set-ups [12-14]. All lengths are rescaled by £q as was explained in Sec. III. Except 
when explicitly mentioned, all results are obtained by using the fully disordered phase of 
the BCP film, (f>A(f) — 0s( r ) — 0.5, as initial condition. Then, a temperature quench is 
performed from the disordered state above the ODT to temperatures below the ODT where 
the lamellar phase is stable. 

A. Chemically striped surface 

The system is modeled using a SCFT scheme for two separate set-ups that are motivated 
by recent experiments [10-14]. In the first set-up, the BCP film is spread on a flat but 
chemically patterned solid surface, while the second bounding surface is the free film/air 
interface, which is either neutral or has a slight preference towards one of the two BCP 
components. In our calculations we take this top surface to be always neutral. A top view 
of the bottom patterned surface can be seen in Fig. 3, and is composed of infinitely long 
stripes in the y-direction of width u s ~ 100 nm that prefer the A component (5u > 0). 
These stripes are separated by neutral inter-stripe regions of width u n having the same 
affinity for A and B (Aw = 0). As the stripes are infinitely long in the y-direction, the 
chemical surface pattern has a one-dimensional square-wave shape and is periodic in the 
x-direction, Au(x + d) — Au(x), with periodicity d = u s + u n 

Au(x) = u s for < x < uj s 

Au(x) — u n for u s < x < d (13) 

Note that we can write formally the surface preference field Aw = u^ — ub as Au(r) = 
Au(x)S(z), where S(z) is the Dirac delta function. All numerical values of Aw are given 
hereafter in terms of its rescaled units, Aw — > NAu. 

In the following we fix the width u s to be twice the natural periodicity, yielding 
w s =2£o=100nm. The phase diagram shown in Fig. 4 is calculated in terms of the film 
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FIG. 3: Top view of a striped surface lying in the x — y plane. The periodicity is d = oj s + Wn, where 
the A-preferring stripes have a width of uj s and the neutral inter-stripe regions are of width ui n . 

thickness, L, and the bottom surface preference, Au, for this set-up and compared with the 
one in Fig. 2 for homogeneous surfaces. All parameters here are taken to be the same as 
for the homogeneous surface, except that the bottom surface has chemical stripes. Further- 
more, we fix the value of the inter-stripe distance (where there is no preferred adsorption), 
to bj n = 8£o = 400 nm so that the pattern periodicity is d = u s + u n = 500 nm, or 10£r> The 
phase diagram is obtained by starting as an initial guess from the perpendicular lamellar 
phase (L_l) or the parallel one (Lm). After convergence, their free energies is compared. 
From the figure it is evident that the Lj_ phase has a larger stability range for the chemically 
striped surface as compared with the homogeneous surface, although the effective value of 
Am on the entire patterned surface is smaller as its value should be averaged over both the 
striped and inter-stripe regions: Au c s = Au{uj s /d). Note that the stability of the Lj_ phase 
is in particular enhanced for special values of L: |£ , f^o 

Due to the existence of many metastable states in BCP melts, the numerical procedure of 
free energy minimization is sensitive to the initial conditions. Instead of converging always 
to the true equilibrium structure at any point of the phase diagram, different metastable 
structures can be obtained. We show some results to illustrate this scenario in Fig. 5. 

For L = 0.7£o, d = 10£o, Au = 1 and Nx = 20 (a typical set of parameters that is located 
inside the L^ stable region), we start with parallel lamellae, perpendicular lamellae and the 
fully disordered state as three different initial conditions and perform a temperature quench 
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FIG. 4: Phase diagram in terms of the film thickness L vs the surface preference An for chemically 
striped surface (dashed line) and homogeneous surface (solid line). The lines separate the parallel phase 
(L||) for larger Ait values from perpendicular one (LjJ at smaller Am values, although the effective Ait 
on the striped surface is Ait e g = Au{ui s /d). The parameters used are Nx = 20, £o = 50 nm and for 
the striped surface: co s = 2£q, uj n = 8^o so that d = 10£o = 500 nm. 

to a temperature below the ODT. The L|| and Lj_ phases in Fig. 5a and 5b, respectively, 
result from quenching from L|| and Lj_ initial conditions. Hence, the system retains its 
orientation after the temperature quench. On the other hand, for the fully disordered initial 
condition, we obtain a mixed structure containing domains of the L|| and Lj_ phases. This 
structure is shown in Fig. 5c and is coined as Lm- 

As explained in Sec. Ill, the maximal deviation of the incompressibility condition, |1 — 
4>A(hJ) ~<t ) B{hJ)\, serves as our accuracy criterion. It is 1.10 x 10~ 6 for the parallel lamellae 
as initial condition (Fig. 5a); 2.57 x 10~ 5 for the perpendicular lamellae as initial condition 
(Fig. 5b); and, 1.07 x 10~ 2 for disordered state as initial condition (Fig. 5c). For the L|| 
and L± it is quite small, yielding a value of about 10 -5 . However, it is not as good in the 
mixed Lm structure (10~ 2 ), because of the existence of internal boundaries between parallel 
and perpendicular domains. To answer the question of metastability we calculate the free- 
energies per chain and obtain /|| = 4.272 > /m = 4.122 > f± = 4.061, corresponding to the 
L||, Lm and Lj_ phases, respectively. Clearly, the most stable structure is the perpendicular 
one, Lj_, and is consistent with our phase diagram in Fig. 4. Note that the free energy 



11 



^ 
N 







0.7 

0.6 

0.5 H 

0.4 

0.3 

0.2 

0.1 

0.0 










0.7 1 

0.6 

0.5 

0.4 H 

0.3 

0.2 

0.H 

0.0 










xlt 



10 15 

xlt 



— I ' 1 — 

10 15 



xlt 



[(b) 



20 



P 



20 



FIG. 5: BCP lamellar structures obtained from numerical solutions of SCFT equations for three different 
initial conditions: (a) parallel lamellar (Lm); (b) perpendicular lamellar (Lj_); and, (c) fully disordered 
state developing into a mixed morphology (Lm)- The film thickness is L = 0.7^o. the top surface is 
taken as neutral, Au = 0, while the bottom one has a striped pattern as in Fig. 3 with Au = 1. The 
inter-stripe widths are set to be co s = 2£q and u n = 8^o, yielding d = uj s + u n = W£q = 500 nm. The 
other parameters are N\ = 20 and £q = 50 nm. The color code corresponds to the four intervals of 
local monomer density < 4>A(f) < 1, as is depicted in part (a). 
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FIG. 6: Calculated BCP lamellar structures for patterned surfaces of increasing inter-stripe distance 
while uj s = 2£q remains fixed (and, hence, increasing d). In (a) uj n = 6£o and d = 84) = 400 nm; in (b) 
u n = 18^o and d = 20^o = 1/im. All lengths are scaled with the lamellar periodicity £q = 50 nm. The 
other parameters are: L = 0.7£o, Au = 1, and N\ = 20. The initial condition is chosen as the fully 
disordered state, ^(r) = 0.5 and all other parameters and color code are the same as in Fig. 5. The 
system exhibits mixed Lm morphologies with Ln regions just on top of the surface stripes and perfect 
L± domains in between the stripes. 

differences between the various states are very small, on the order of 2-5%, manifesting the 
tendency of the system to get trapped into metastable states. 

Our findings have also experimental implications because in experiments the film struc- 
ture depends strongly on its history and sample preparation [14, 33, 34]. The claim is that 
once the system is prepared in its L^ it will stay there. But if the film is prepared above 
the ODT, in its fully disordered state, the film can get stuck in a metastable mixed lamellar 
structure, Lm- Although in experiments it is not always possible to heat the system above 
its ODT because of polymer break-down and oxidation, in many cases, higher temperatures 
are used to anneal the film and allow it to reach its final state via faster dynamics. 

Another interesting feature is presented in Fig. 5. Perfect perpendicular lamellar struc- 
tures between neighboring stripes are visible. Furthermore, we can obtain such perfect L^ 
structures for a wide range of small and large periodicities, ranging from d = 400 nm in 
Fig. 6(a) to d — ljum = 20£q in Fig. 6(b). However, we find that it is difficult to get rid 
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of the parallel lamellar regions induced by the striped pattern, even when we further re- 
duce the BCP film thickness L to values much less than £q and decrease the values of Au. 
Furthermore, a preliminary study [49] indicates that slow temperature annealing from the 
disorder state (above ODT) to the ordered lamellar state (below ODT) does not seem to 
prevent the formation of the mixed Lm phase. 

B. Periodic grooved surfaces 

In order to overcome the problem of getting trapped in Lm mixed states and inspired 
by recent Nanolmprint lithography (NIL) experiments [14, 32-34], we explored yet another 
type of surfaces. The set-up can be seen in Fig. 7, where the BCP film is confined between 
two solid surfaces. The bottom surface at z = is flat, while the top one has a periodic 
arrangement of grooves (along the x-direction) made of a series of down-pointing 'fingers' 
of thickness ui, separated by inter-grooves regions ('plateaus') of thickness Uh- The periodic 
height profile h(x) = h(x + d) has the form: 

h(x) = L[ for < x < u>i 

h(x) = L h for ui < x < d (14) 

where the height is measured from the z = surface. Formally, Au(r) = u& — ub used in 
the solution of Eqs. (6)-(7) is given by Au(r) = Au(x)S(z — h(x)). 

The figure shows the surface height profile in the x — z plane, for profiles that are trans- 
lationally invariant in the y-direction. The periodicity in the x-direction is d = ui + Uh 
and the finger width is chosen to be ui = 5£q = 250 nm. The top surface (mold) is put 
in direct contact with a BCP film spread on a neutral and flat bottom surface (at z — 0). 
The distance of closest approach between the two surfaces is Li, while the maximal height 
difference between them is Lh- This means that the finger height of the mold is Lh — L\. 
Assuming film incompressibility, we get a relation between the thickness L of the original 
BCP film and the two height parameters, L\ and Lh- L = (LiUi + Lh^h)/d. In experiments, 
the average thickness L is fixed, while in the numerical study, we control directly L\ and Lh- 

By varying the values of the parameters d, Lh and L\ of the mold, and the strength of 
surface interactions Au, we can get a sequence of BCP patterns. Furthermore, we obtain 
perfect perpendicular lamellar structures extending throughout the film thickness for some 
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FIG. 7: A cut (side view) through the top grooved surface (the mold). The periodicity in the ^-direction 
is d = u>i + uJh, with ui and u>h being the finger and inter-finger width, respectively. L\ is the distance 
of closest approach to the bottom surface at z = 0, and Lh is the largest film thickness. The initial 
film thickness is equal to the average film thickness in the mold, L = (L[Ui + LhtOh)/d. 

special patterned surfaces. 

We calculate the phase diagram in terms of the maximal film thickness, Lh, vs the surface 
preference, Am = ua — ub- The interaction strength on all exposed surfaces of the upper 
grooved mold have the same value of Au. In addition, we set Li = 0.3^o, d = 10^o and 
ui = 5^o- The result is shown in Fig. 8, from which we can infer that this set-up greatly 
affects the phase diagram as compared with Fig. 2 for a uniform Au surface. The transition 
line from L|| to L^ is shifted upwards so that its minimum is obtained for Lh = £ Q where 
Am = 4.4. This is similar but more pronounced than the behavior seen in Fig. 4 for the 
chemical striped surface around the L/£ = 1.0 region. 

However, when we start from a fully disordered state as initial condition inside the stable 
L± region of Fig. 8 (e.g., Lh = 0.8£o and Am = 0.1), we do not get the fully perpendicular 
lamellae L^ but rather a mixture of parallel and perpendicular lamellar regions (the Lm 
structure) as shown in Fig. 9. 

We find two ways to improve on the perpendicular orientation by changing the mold 
geometry and surface characteristics. First, we decrease the film thickness by decreasing Lh 
to 0.6^0- In this case, we do a gradual temperature quench, starting from the disordered 
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FIG. 8: Phase diagram in terms of the maximal film thickness L^ vs the surface field difference of the 
two blocks, An, for a periodic grooved surface (dashed line) and a homogeneous surface (solid line). In 
the latter case, the film thickness L is equated with Ly t . Other parameters are: £q = 50 nm, L\ = 0.34, 
d = 154 = 750 nm, ui = 54 = 250 nm, co h = 104 = 500 nm and N X = 20. 
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FIG. 9: BCP density distribution for d = 154 = 750 nm and N\ = 20 starting from a fully disordered 
initial condition. The bottom surface is neutral and the top surface has a square wave height profile as 
in Fig. 8 where An = 0.1. Other parameters are L^ = 0.84 and L\ = 0.34, ^h = 104, wi = 54, 
yielding L/£ = 1.9/3 ~ 0.64. 
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FIG. 10: BCP density distribution for d = 15£o = 750 nm and N\ = H-5 when the initial condition 
is the fully disordered state. The system is first annealed to N\ = 11-5 in (a) and then to N\ = 20 
in (b). The bottom surface (z = 0) is neutral and the top surface has a square grooved structure with 
Au = 0.1. Other parameters are Lh = 0.6^o, Li = 0.3^o, ^h = 10^o ar| d ui = 5^o, yielding L = 0.5^o- 



state above the ODT and quenching to temperatures just below the ODT, N\ = 11.5, 
and only then proceed with a deep quench to N\ = 20. This two-step procedure is shown 
in Fig. 10(a) and (b). Perfect perpendicular lamellar structures emerge. Moreover, using 
this two-step procedure we can even obtain a perfect perpendicular lamellar structures with 
much wider u n yielding d = 1.25 /im (or equivalently d/£o = 25), as is shown in Fig. 11. 

A second variation is to construct the grooves from two separate materials with different 
A/B preference. The protruding 'finger' parts are assumed to have a small A preference 
(Am = u\ > 0) both on their vertical and horizontal parts, while the high plateau parts are 
taken as neutral (Au = 0). 



Au(x) = ui for < x < u s 
Au(x) = for ou« < x < d 



(15) 



With this special surface geometry and interactions, we obtain perfect L^ structures for 
wide range of film thicknesses. An example for such a set-up with d = 25£q periodicity is 
shown in Fig. 12. 
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FIG. 11: BCP density distribution for d = 25^o = 1-25 fim using two-step annealing procedure. First 
to N\ = 11-5 in (a) and then to N\ = 20 in (b). The initial condition is the fully disordered state. 
Other parameters are: L^ = 0.6^o, <^i = 5^o. ^h = 20^o and L/ = 0.3£rj, yielding L = 0.54£o- 

V. DISCUSSION AND CONCLUSIONS 

In this paper we addressed several surface patterns as inspired from recent experiments 
in relation with ordering and orientation of lamellar phases of block copolymer (BCP) films. 
In the first set-up, we model a BCP film confined between a chemical striped solid surface 
and the free film/air surface. In a second set-up, the film is considered to occupy the gap 
between two solid surfaces; a flat one and a hard mold with specific square-shape grooves as 
is inspired from recent Nanolmprint lithography (NIL) experiments. 

The main question both experiments and modeling should attempt answering is how to 
induce a perfect perpendicular order in BCP films? In particular, how this can be achieved 
using patterned surfaces with structural features (stripes and grooves) that have a periodicity 
d much larger than the lamellar periodicity £q. Having such sparse surface features will 
reduce substantially the cost of large-scale production of surface templates and BCP films 
and is essential for applications, e.g., in microelectronic and nano- lithography processes. 

Using the first set-up of the chemically stripes on an otherwise flat and neutral surface, we 
are able to show that the perpendicular phase L^ has a larger stability region in parameter 
space described by the film thickness L and surface preference (A-u), as compared with the 
homogeneous surface. Note that this is the case even for inter-stripe distances u n that are 
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FIG. 12: BCP density distribution for d = 254 = 1.25 /zm. The bottom surface at z = is neutral, 
Au = 0. The top surface is a square grooved with Au = 0.1 on the sides and tips of the grooves and 
neutral (Au = 0) on the top plateau parts (see text) and Eq. (15). Other parameters are Nx = 20, 
L h = £q, Li = 0.34, oj h = 204 and u>i = 54, yielding L ~ 0.864. 



an order of magnitude larger than the stripe thickness, u s . This is in spite the fact that 
the effective (averaged) Aw for the striped surface is smaller than the corresponding Au 
on the homogeneous surface, Au c g = Au(u s /d) < Au. We equally find that the system is 
very sensitive to initial conditions. Starting from a fully disordered state, above the order- 
disorder temperature (ODT) and annealing the temperature into the lamellar region, will 
mainly produce a mixed morphology Lm as can be seen in Fig. 5(c) and Fig. 6. Although 
the stripes nucleate growth of BCP layers on top of them (namely, domains with parallel 
orientation, L||), perfectly oriented perpendicular domains, Lj_, are induced on top of the 
neutral inter-stripe region. 

In our model the Lm mixed morphology is a result of the large number of metastable 
states (local minima) which the system possesses. Although the true equilibrium is the L^ 
phase, it is hard to find it numerically unless one starts with the proper initial conditions. 
This drawback should also be expected in experiments, where during sample preparation, 
the film undergoes many external stresses and defects are abundant. It will be of interest 
to verify in experiment our findings by doing a slow temperature annealing of BCP films 
from their disordered liquid state (above ODT) into the lamellar region (below the ODT). 
Such a slow temperature annealing has the potential to produce highly oriented BCP films. 
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Although not in all systems it is possible to reach temperatures above the ODT without 
damaging the BCP chains, we equally note that in many cases annealing at high enough 
temperatures has the advantage that the system can reach it final state with faster kinetics. 

In the second set-up, we modeled a hard mold which is pressed onto a BCP lamellar film. 
We show that this Nanolmprint lithography (NIL) process greatly enhances perpendicular 
order in lamellar phases. Perfect Lj_ can be seen for film thicknesses below £$ even when 
the groove width Uh (filled with the BCP film) is five times (or even larger) than the solid 
'finger' (ui) sections. Here the slow annealing from above ODT to below the ODT is very 
successful, demonstrating that this set-up is more suitable for lamellar orientation purposes 
than the chemical stripe set-ups discussed above. 

In Fig. 12 we proposed a mold with even superior orientation qualities. For this mold the 
surface preference of the downward protrusion sections (the 'fingers') is larger than that of 
the top section of the groove (plateau-like). As the latter preference interferes with the L|| 
ordering, reducing this surface preference will enhance Lj_ ordering, especially in the desired 
case of thin fingers and wide plateaus, where Uh ^> 0Ji- 

In experiments, it is harder to produce a mold with such specific surface characteristics as 
seen in Fig. 12. One way would be to form it from two separate materials or to use a selective 
coating during mold preparation. However, creating such a mold can be a costly and delicate 
process that will be hard to mass reproduce. Yet another possibility is to have an effective 
chemically heterogeneous mold shown in Fig. 13. Suppose that the groove height is only 
partially filled with the BCP melt, creating pockets of air on the top of each groove [32]. 
The film/air interface within each groove can be thought of as another interface with almost 
neutral preference the two blocks. This situation amounts to taking different values of Au on 
the finger-section and plateau-section of the mold [see Eq. (15)]. While Am = on the top 
section (plateau) of the groove, it is non-zero on the mold 'finger' sections. This is exactly 
the situation explored in our calculations and shown in Fig. 12, and it may be worthwhile 
to further explore this partial filled mold in future experiments. 

Our theoretical modeling relies on numerical solutions of self-consistent field theory 
(SCFT) equations. We minimize the corresponding free energies and converge to film mor- 
phology whose free energy is an extremum using an iterative procedure. We find that the 
numerical procedure is sensitive to what is used as initial conditions for the BCP structure. 
The convergence can be towards local (metastable) states and not always towards the true 
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FIG. 13: Schematic drawing of the Nanolmprint lithography (NIL) set-up where the mold only partially 
is filled with the BCP film. Effectively, this means that the film sections in contact with the side 
boundaries of mold feel a different surface field than the top horizontal facets, which are exposed to 
the air. 

equilibrium. This is an unavoidable feature of the numerical procedure. It is not an artifact 
but rather reflects the true physical situation as seen in experiment. The BCP film has many 
metastable states separated by energy barriers and it is hard to reach the true thermody- 
namical equilibrium state. Slow annealing from above the ODT or from high temperatures 
is one way to overcome this difficulty, at least in a partial way. 

It will be of great interest to further proceed and extend our two-dimensional calculations 
to full three-dimensional ones. This will require much longer computation times but will 
allow us to distinguish between perfectly oriented perpendicular lamellae and those that 
stand up but which also wander around in the x — y plane. For applications it is important 
to have perfectly oriented Lj_ phases in the z-direction, that are well aligned in the lateral 
(in-plane) directions. 

Although our present study is not exhaustive, it shows many possibilities of explaining 
some of the experimental findings and even points towards interesting directions for future 
experiments. 
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